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Abstract 

^SJ ' We investigate the possibility that the large mixing of neutrinos is 

^ . induced by their large coupling to a five-dimensional bulk neutrino. In 

the strong coupling limit the model is exactly soluble. It gives rise to 
an oscillation amplitude whose squared-mass difference is independent 
. of the channel, thus making it impossible to explain both the solar and 

I the atmospheric neutrino oscillations simultaneously. 

O ; 

PhI 1 Introduction 

O . The idea of using an extra dimension to explain physics goes back to Kaluza 

and Klein |^. Spurred by the discovery of D-branes in string theory [Q], 
this idea was revived in recent years by introducing branes P] with many 
^ ■ variations to the theme. We will adopt here the simplest one, where a sin- 

■ gle three-brane is embedded in an extra fiat dimension. All Standard-Model 

(SM) particles are confined to our 3-1-1 dimensional world, the three-brane. 
Gravity, right-handed neutrinos, and other particles without SM quantum 
numbers are allowed to live in the 4-1-1 dimensional bulk. This modification 
thus introduced to gravity has been tested experimentally [|]. Its conse- 
quence on astrophysics and cosmology is extensively investigated In 
this paper we will discuss an application of this picture to neutrino physics 

If solar and atmospheric neutrino deficits are attributed to neutrino oscil- 
lations, then neutrinos must be massive. Direct measurement puts an upper 
bound of about 2 eV on [^. If it is a Majorana particle, the failure to 
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detect neutrinoless double j3 decay further limits a combination of the light- 



est neutrinos and mixings to have a mass less than 0.26 eV |11]. In any case 



neutrino mass is much smaller than other fermion masses and this mystery 
should be understood. Moreover, the amount of mixing needed to explain the 



atmospheric neutrino deficit is large ]12[, whereas the mixings of quarks are 



small. Presently these two differences between neutrinos and other fermions 
are not clearly understood. 

In the seesaw mechanism, the smallness of neutrino mass is explained by 
the high energy scale of the right-handed neutrinos. The smallness can also 
be attributed to the large size of the extra dimension f^, 0, |^ , and it would 
be nice if the extra dimension can also explain the large neutrino mixing. It 
is this latter possibility that we carefully examine in the present paper. 

To this end we study the simplest model, which consists of one bulk 
neutrino and three active Majorana brane neutrinos, z/g, z/^, and Vr- Since 
mixings between quarks are small, to simplify matters we shall assume no di- 
rect mixings between the active neutrinos in this model. All mixings between 
them are induced solely through their mass couplings to the bulk neutrino. 
The bulk neutrino is a Dirac fermion in 4+1 dimensions, and a SM singlet. 
All in all, the theory contains six parameters, the Majorana masses nii for 
the three active neutrinos, and their mass couplings di to the bulk neutrino. 
For equal couplings this model was discussed by Dienes et al in 0. [Other 
similar models were also studied [0].] Without Majorana masses and some- 
times with only one generation, this model has also been studied by several 
other groups 0. The mass couplings di modify neutrino masses and induce 
a mixing between them. If the large mixings seen in experiments are to be 
explained by this simple model, the mass coupling strengths must be large. 
We recall that di = UiV / yj2'KM^R where is the Yukawa coupling between 
the i-th brane neutrino and the bulk neutrino, v is the usual vacuum ex- 
pectation value of the SM Higgs field, is the fundamental 5D scale and 
R is the compactification radius. Hence, the strong coupling case will arise 
from large brane- bulk coupling and usual perturbative treatments assuming 
di ^1 will not apply. It is this strong-coupling scenario which we find most 
interesting and will be the main focus of this paper. 

Numerical studies have been carried out in |^, 0] for specific parameters. 
Our analysis is completely analytical and we find remarkably that there is 
actually an exact solution to the problem in the strong-coupling limit that we 
are interested in. In that limit, we are able to give a complete discussion of 
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the eigenvalues and eigenvectors, as well as exact brane-to-brane transition 
amplitudes which we have not found in the literature. 

In the absence of coupling, the mass eigenvalues are integers for the in- 
finite number of sterile neutrinos composing the bulk neutrino, and rrii for 
the active neutrinos confined to the brane. All masses are understood to be 
measured in units of the inverse size of the extra dimension. In the strong- 
coupling limit, sterile and active components are thoroughly mixed up. The 
mass eigenvalues now consist of every half-integer, plus two others, Xb and 
Ac, determined by and the fixed ratios of the rfj's. More generally, if there 
are / active neutrinos then there are / — 1 such other eigenvalues. With 
this knowledge, the transition amplitudes from active flavor i to active flavor 
j = z/e, z/^, z/,-) can be computed exactly at any 'time' r. It contains 
three terms. The first comes from a sum of all transitions via half-integral 
eigenvalues. It is short lived, and it decays to zero within a short 'time' 
inversely proportional to the fourth power of the coupling. The second and 
the third terms are transitions mediated by eigenvalues A;, and Ac. After the 
transient period when the first term dies away, we are left with an amplitude 
resembling the usual two-neutrino oscillating amplitude, but with a squared- 
mass difference Ac — A^ independent of i and j. This independence makes 
it impossible to explain both the solar neutrino deficit {i = j = Ve) and the 
atmospheric neutrino deficit {i = j = i^p) simultaneously, for their measured 
squared- mass differences are at least two orders of magnitudes apart. There 
are actually also other reasons why this model would not work in practice 
but they will be discussed later in the main text. 

The failure of this model to describe Nature is not due to the presence of 
sterile neutrinos per se, for they are still allowed. Although sterile neutrinos 



are disfavored in the solar and atomospheric data of Super-Kamiokande [12 



solutions still exist in which the solar oscillation occurs between z/g and a 
sterile neutrino 0, and a substantial sterile component is also allowed in 
atmospheric neutrino oscillations [O]. If LSND experiment is confirmed we 



do need a sterile component anyway. The failure rather is due to the presence 
of an infinite number of them, with mass eigenvalues at half integers. It 
is the regular spacing of these eigenvalues at large coupling that produces 
a destructive interference at large r, which eliminates the first term of the 
transition amplitude, and kills the model. In this our conclusion is consistent 
with earlier numerical investigations and analytical studies for one active 
neutrino P, HI- If the first term were oscillatory instead of decaying, which 
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would be the case if the infinite number of them were replaced by a single 
sterile neutrino, then we would just get a regular three-neutrino oscillation 
formula and a problem of this nature would not arise. 



2 Mass Matrix and Its Diagonalization 

We consider a model in which a three-brane is placed at an / Z2 orbifold 
fixed point. The size of the extra dimension is measured by the radius R of 
the circle S^. The model consists of three left-handed Majorana neutrinos 
and z/t, confined to the brane, and a Dirac bulk neutrino ^ living in 
the 4+1 dimensional world. The only interactions are mass terms, and we 
shall express all masses here in units oil/ R. The bulk neutrino decomposes 
into an infinite number of 3+1 dimensional sterile neutrinos, with masses 
77, G Z. The Majorana masses of the brane neutrinos are denoted by mj, and 
their couplings to the bulk neutrino are denoted by the complex number di. 
Since phases are not currently detectable, we will assume rrii and di to be 
real. This practice is also adapted in numerical studies 0, |^. 
The neutrino mass matrix is given by H] 
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Let V = (fi, V2, ■ ■ ■)'^ be its column eigenvector with eigenvalue A. Let A = 
Y1^4: "^j' ^"^^ ^ — Si=i diVi- The eigenvalue equation M-v = Xv in component 
form is 



TUiVi + diA 
b + CiVi 



Xvi 



(z = 1,2,3) 
(^>4) 



(2) 
(3) 
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where C2n = 2 — n and C2n+i = n — 1 for n > 2. It follows from (Q) that 

b 

b 

^2n+3 — "T"; ! 

X + n 

V2n+4 = ^ ^ ^ , > !)• (4) 



Hence 



m[ + 


did2 


did3 \ 


u 


d2di 


m'2 + d\ 




■ V2 


d^di 


dzd2 


m'^ + dl) 


\V3 



A = Vt;, = 6A V — ^ = (5) 
^ ^ - tanfvrA) 

i=4 n=— 00 ^ ' 

Substituting this into (|), it becomes 

did-^ \ I vi\ 

= 0, (6) 

\ d^di d3d2 m'^ + '^3 / J 
where 

= (mj — A) tan(7rA)/7r. (7) 

The eigenvalue A is then determined by det[M'] = 0. We shall use the 
characteristic equation in the form 

i=l ' 



2.1 Eigenvalues and Eigenvectors 

The eigenvalues can be determined in the following way. In the unit interval 
(n— |, n+ 1), n G Z, the left-hand side of increases from —00 to +00 as A 
moves from the left edge to the right edge of this interval. On the other hand, 
the function on the right-hand side of is a monotonically decreasing func- 
tion of A in each of the four intervals (—00, mi), (mi, 7712), (m2, ms), (ma, +00), 
if mi < m2 < m3. It decreases from to —00 in the first interval, from +00 
to —00 in the second and third intervals, and from -|-cxd to in the fourth 
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interval. Hence if we divide the real axis into intervals bounded by half in- 
tegers (odd integers divided by 2), as well as mi,m2,m3, then there is one 
and only one solution of (|^) in each of these intervals. See Fig. 1. 

From (H), we can determine (f i, f2, fa) simply by taking the cross product 
of any two rows of the matrix M'. Using (^ and (|^), the result (to within 
an arbitary normalization) can be taken to be 

Vi = -^, (^ = 1,2,3). (9) 



With this normalization, 

b = Y^ diV, = — = - tan(7rA) , (10) 



X — rrii 
i=i 1=1 



where (H) is used in the last expression. The rest of the components, Vi for 
i > 4, is determined by (^. The norm of the vector v is equal to 



di ^ ^ 



^\X-miJ cos2(7rA) ^ ^ 



The second term in (|TT|) is obtained from 



i=A 



(A^ — n^)^ \sin(7rA) 

by using (|TUp. 



2.2 Unitary matrix 

It is more convenient to replace the flavor indices i = 1,2,3 by 01,02,03. 
Index i = 2n + 3 for n > 1 will be renamed index —n, and index i = 2n + 4 
for n > will be renamed index n. With this alteration, eqs. (|^) and (Q) 
become 

(01,02,03) = i'^c^f^'^r) 

{n G Z). (13) 



_ di 
X — rrii 
b 

X — n 
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We shall refer to Oj as the brane components, and n as the bulk components. 

The overlapping matrix between a normalized eigenvector |A) and a flavor 
component (/| (/ = ai,n) is then given by the unitary matrix 

Ufx = (/|A) = ||^, (14) 
where f/(A) is given by ([13|) and -B(A) is obtained from (|ll|). 



3 Neutrino Oscillations 

The probability for a brane component aj to turn into a brane component 
Oj, after the neutrino of energy E has travelled a distance L, 

V^,{r) = \A,ir)\\ (15) 
is given in terms of the transition amplitude 

A,(r) = J2KxUa,xe-'''\ (16) 

A 

where r = L/2E for vacuum oscillations. The sum over eigenvalues A is an 
infinite sum and can be evaluated only when the eigenvalues are explicitly 
known. 



4 Strong Coupling Limit 

When di = 0, the neutrino mass matrix M is diagonal. The eigenvalues are 
rrii for the brane states of active neutrinos, and G Z for the bulk states of 
sterile neutrinos. There is no mixing between any two flavor states. 

When di ^ 0, eigenvalues shift, and mixing between flavor states are 
induced. If di are small, the shifts and mixings can be calculated by per- 
turbation theory. This however is not the interesting regime because it will 
never generate the observed neutrino flavor mixings, which are large. 

When di are big, there is a chance that such large mixings can be induced. 
This is the regime we will study in the present section. Remarkably, in this 
strong coupling regime, the eigenvalues can be obtained explicitly and the 
infinite sum encountered in (p!6| ) can be computed. 
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4.1 Eigenvalues 

Let di = dci, J2'i=i^l — 1; d ^ 1. For fixed Cj and increasing d, the right- 
hand side of (H) increases in magnitude, so the left-hand side must change 
accordingly. If the equality in (|^) is positive, the eigenvalue A must move to 
the right (R). If it is negative, A must move to the left (L). 

Divide the real axis into six sections, Ji = (— oo, mi), J2 = {mi, \b), J3 = 
{Xb,m2), J4 = (m2,Ac), J5 = (Ac,"^3), Je = {^,3,00), where A;, and Ac are 
the two zeros of 

3 2 



r(A) = $^-^, (17) 



=1 



explicitly given by 

p = e^(m2 + ms) + 62(7713 + mi) + 63(7711 + m2), 

c = [p^ — 4(6^m2m3 + 62m3mi + 63mim2)] ^ . (18) 

Then A moves L if it is in Ji,3,5, and A moves R if it is in J2,4,6- See Fig. 1. 

Remember that if the real axis is divided into intervals I\ bounded by 
the half integers and m^, then there is one and only one eigenvalue A in each 
interval. If I\ C ^1,3,5, then A moves L. If Ix C ^2,4,6, then A moves R. If 
Ix contains either Af, or Ac, then A moves towards these zeros of 7'(A). In 
particular, if m^ is one of the two boundaries of Ix, then A always moves 
away from that boundary. If both boundaries are some m^, then a zero Ab,, 
must be contained in that interval where the eigenvalue moves towards. In 
any case, the eigenvalue always stays away from the m^ boundaries. 

As d 00, the eigenvalues will therefore end up either at the half- integer 
boundaries, or a zero A^, or Ac- In fact, there is one and only one eigenvalue 
at every half integer. 

To see that, choose any half integer, and denote the interval to its left 
(right) by /_ (/+). The only way this half integer ends up not to be an 
eigenvalue, or end up to be equal to two different eigenvalues, is when the 
two A's in J_ and /+ move in opposite directions as d increases. This however 
is not possible unless a zero is contained in /_ or in /_|_. If it is in J_, then 
A G /- will drift towards this zero, and A G /+ will move L, drifting towards 



c 
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the half-integer value in question. Similar thing happens when the zero is in 
J+. We therefore conclude that there is one and only eigenvalue at each half 
integer value. 

It is important to note that the eigenvalues A;, and Ac are special simply 
because they do not fall into this infinite sequence of half integers. It is not 
because they are brane states. In the strong coupling limit, brane and bulk 
components are thoroughly mixed up. Nevertheless, we may still define a 
brane state as one that extrapolates back to the eigenvalues when di = 0, 
and a bulk state as one that extrapolates back to an integer eigenvalue in 
that limit. With that definition, each zero may correspond to either a brane 
state or a bulk state. If Ix containing this zero also contains an integer, then 
it is a bulk state. If it does not, then it must be a brane state. 



4.2 Transition amplitude 



With this strong-coupling spectrum, the transition amplitude in (^6[) can be 
written more explicitly as 



n=— oo 



+ 



+ 



5^(A„)(A„ - mi)(A„ - rrij) 



B^{\){h - "^i)(A6 - rrij) 
B^{Xc){Xc - rni){Xc - rrij) 



(19) 



For (i ^ 1, there is a qualitative difference between the normalization 
factors -B^(An) and i?^(Afe,c)- This is so because the right-hand side of (H), 
which according to (|17D can be written as (i^r(A), behaves differently for the 
two sets of eigenvalues. For A^^c, ''"(A) = Oil/dP) so these eigenvalues as well 
as cos(7rA) are finite in the limit d oo. For A„, r(A) = 0{d'^/{X — rrii)), so 
in the strong coupling limit tan(7rA) as well as 1/cos{7tX) are of order c?^ for 
finite A, and are of order 0{d'^/X) when A is comparable or larger than d"^. 



Using (11), it follows then that B^{X) is of order d^ for Xb.c, but is of order 
d^/ X for An. 

Returning to (|l^), it follows from this observation that the second and 
third terms on the right are finite in the strong coupling limit, but each of 
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the terms in the infinite sum is of order l/d? when A„ <^ (P. This means 
that we can freely drop a finite number (say N) of terms from the infinite 
sum in the strong couphng hmit, and retains only those terms in the infinite 
sum with large A„. In that case, the right-hand side of (||) becomes d?/\, 
and we can use it to solve for 1/ cos^(7rA) to get -B^(A) ^1/ cos^(7rA) = 
1 + + Ti^(f)/\^ = 1 + K^/\^. Since AA„ < (P and N/(P < 1, we may 
now replace the infinite sum on the right of (|l^) by the integral 

At r = 0, the integral can be evaluated to give /(O) = Tr/i^ ~ l/dP. 

Moreover, for ^ 1, it follows from (0) that B^{X) ~ (Ps{X) is true for 
A = Afe, Ac (assuming that cos(Ab^c) 7^ 0), where 



Therefore, 

Aij{0) = CiCj 1 + 



1 



s{Xb){Xb - m){Xb - rrij) 
1 



+ 



s(Ac)(Ac - mi) (Ac - mj)_ 



(22) 



According to (|T^), this should be equal to 6ij for all parameters Cj and rrii, 
if the infinite-dimensional matrix Uf\ is unitary, as expected from the com- 
pleteness relation. This is indeed so and can be verified numerically. To 
prove the identity analytically, a change of parameters is desirable. 

The model is specified by six parameters: di = dci and rrii. They are 
actually not the most convenient parameters to use because according to (|^), 
a shift of the m^'s by a common integer simply shifts all the eigenvalues by 
the same integer, so the more useful parameters are X — rrii rather than rrii. 
In fact, all the quantities in ( p2D can be expressed in terms of the parameters 

1 1 

Xi = T , Vi = T . (23) 

Xb -rrii Ac - rrii 
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Note that these two sets of parameters are related by 

- - - = Ae - Afe = c, (24) 

SO that there are only four independent parameters in the set: c and Xj, or c 
and Hi. 

We can use the three conditions, = 1 and r{\b) = r(Ac) = 0, to solve 
for the three variables , Cg, e|. The result is 

e? = {x2y3-y2X3)/D, 
4 = {x^yi - y^xi) / D , 
4 = ixiy2-yiX2)/D, 

D = {xiy2 + X2y3 + xsyi) - {yiX2 + y2X3 + y3Xi). (25) 
We can use (|^) to eliminate either yi or Xi, and we get respectively 

Sb = s{Xb) = CX1X2X3, (26) 



or 



— t\ \ X1X2X3 , . 

Sc = s{Xc) = -cyiy2y3 = -cy— r-— TT— r. (27) 

(1 + CXi){l + CX2){1 + CX3) 



In these notations, we can write ([2^ ) in a more transparent form as 
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where 



A,(0) = J2KV^p^ (28) 



ei eixi/y/sl eiyi/^c 
V = \ 62 62X2! (i^y^l^/s'c I , (29) 
.63 63X3/^55; 632/3/^ 



up to an arbitrary phase for each column. The matrix V is unitary, as can 
be verified from the explicit equations (12^) to (|27|) . 

The matrix V is parametrized by four real parameters, c, Xj, plus three 
arbitrary phases of Cj, more than the usual 3x3 mixing matrix. Nevertheless, 
there are mixing matrices not of the form (pQ]), with the constraint (E3) to 



11 



(p7[). An example is the bimaximal neutrino mixing matrix (53 
C3 = cos6'3): 



u 



C3 
-1(1 
1(1- 



1^2 



c^l \f2 S3 

-(I-S3) C3/V2 



HS3j 

S3) 



sin 513, 



(30) 



We are now ready to discuss the case r > 0. The function /(r) can be 
expressed in terms of error function with an imaginary argument, but it is 
more straight forward to write it directly in the form /(r) = {71/ K)g{K'^T), 
where 



9[x) = 



du- 



(31) 



The absolute magnitude of the function g{x) decreases monotonically from 
1 to zero. For large x, g{x) ~ (1 — i)/ \/2'kx. 

The transition amplitude (19), in the strong couphng limit d ^ 1, is 
therefore equal to 



g{K'r) + 



+ 



s(Afc)(Af, - mj)(Ab - rrij) 
s(Ac)(Ac - mi){\c - rrij) 



(32) 



with the matrix V given in (^). At any 'time' r, the amplitude is described 
by six parameters: K, c, Xi, X2, X3, and one of the mj's so that A^ and Ac can 
be determined from Xi and yi respectively. Alternatively, we can let A^ be 
the sixth parameter, whence Ac = A;, + c. 



4.3 Comparison with direct mixing 

In this theory, neutrino mixing is induced by coupling of active neutrinos to 
bulk neutrinos. No direct mixing of active neutrinos takes place. Neverthe- 
less, in the strong coupling limit, the resulting amplitude ( |32| ) bears a close 
resemblance to the usual theory in which sterile neutrino is absent, and mix- 
ing of active neutrinos is directly introduced through the 3x3 MNS matrix 
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u [|I4|. If /Ui,/i2,/i3 are the mass eigenstates, the transition amphtude with 
direct mixing is given by 

= u'^^Ujie-'^''^ + u*^Uj2e-"''^ + u%uj3e-'^'^^. (33) 

Since both V and u are unitary matrices, if we identify with /x^ 
with /i| (modulo a common constant), then the main difference between (^) 
and (^) is the replacement of g^K'^r) by e~*'^i'^. Both functions equal to 
unity at r = 0, but the latter oscillates and the former decays to zero within 
a short transient 'time' ~ l/K"^ <^ 1. The oscillations mediated by three 
mass eigenstates in the direct-mixing theory, is replaced after the transient 
period by an oscillation mediated by two eigenstates (A;, and Ac) in the bulk 
theory. Note that these two eigenstates can be either brane states or bulk 
states (see the discussion at the end of Sec. IIIA), so in general they have no 
direct connection with the three mass eigenstates in the direct-mixing theory. 
In particular, we should not think that one of these three mass eigenstates 
gets lost after the transient period. 

Moreover, unlike the direct-mixing theory in which the total probability 
to be in some active neutrino state is always 1, in the bulk theory this total 
probability decreases with r because of leakages into the bulk. Quantita- 
tively, 

V,{T)^Y.V,,{T) = l-e]{l-\g{K'T)\'}^l-e] (r » i^-^). (34) 

i=l 

As a check we go to the limit of equal couplings, i.e. all e^'s are equal. 
Since Ylii = 1 then each ef = 1/3 and the above probabihty approached 
2/3 in agreement with the result of [^. 

In the bulk neutrino theory, with an infinite number of sterile neutrinos, 
is very different from the direct-mixing theory with three active neutrinos, 
or three active neutrino plus one sterile neutrino. In particular, there is no 
way the bulk neutrino theory can explain both the solar and the atmospheric 
neutrino deficits simultaneously, at least not in the strong-coupling limit after 
the transient period is passed. 

To see that, we use ( P^D to write the survival probability after the tran- 
sient period in the form 

VuM = \AuiT)\' =p,~p[ cos' (^A\ 



13 



p. = {\Vi2\^ + \V.3\Y = {l-ei)\ 
P- = 4(l-e^)|\/,3p, 
= Xl-Xl 



(35) 



This is to be compared with the usual two-term survival probability formula 
for solar and atmospheric neutrinos, 



where for solar neutrino oscillations is at least two orders of magnitude 
smaller than A^ for atmospheric neutrino oscillations. There is no way the 
formula (pSj), with A independent of i, can fit both the solar and the atmo- 
spheric deficits. 

We might want to be less ambitious and use (^) to fit only solar (or 
atmospheric, but not both) neutrino deficit by adjusting the parameters for 
( |35D to emulate (|36| ) with the experimentally measured 6i and Aj = A. Even 
so it will not work, not exactly anyway. The reason is that to make pi = 1, 
we must put Cj = 0. In that case ( [T7| ) has only one but not two zeros, and the 
three terms appearing in (|32D become two terms. After the transient period, 
we are left with one term, whose probability does not oscillate at all. 

These failures are due to the disappearance of the first term beyond the 
transient period. What about taking the opposite limit when r is so small 
that g{K'^T) ~ 1? That would not work either for two reasons. First, in 
order for oscillations to occur we must have A^r and/or A^r to be of order 1, 
which requires A^ ^ and/or A^ ^ K^. This cannot be achieved within 
the strong coupling limit when and rrii are fixed and d ^ 1. Moreover, 
even if we ignore this difficulty, there is no way the unitary matrix V can 
simulate a realistic MNS matrix u, as pointed out previously in (pO|). 

In conclusion, the idea of using a bulk neutrino living in an extra dimen- 
sion to explain neutrino oscillation is a very attractive one. It might explain 
why neutrino masses are small and in the strong brane-bulk coupling limit, 
why their mixings are big. Unfortunately, such strong coupling also causes 
a large leakage into the sterile modes in the bulk, so that this theory in the 
strong-coupling limit is not compatible with experiments. Whether an inter- 
mediate coupling is viable or not depends on whether rapid variations as a 
function of energy and distance are observed in future experiments. 




(36) 
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Our study complements the mostly numerical work done on this simple 
model. Since the solution is nonperturbative, we deem this work to be an 
important first in the study of the physics of strongly coupled bulk neutrinos. 
As more data in solar and atmospheric neutrinos become available there will 
be more serious attempts at seeing whether bulk neutrinos can explain all or 
part of such improved data sets. Our work can be used as a guide to such 
attempts. 

This research is supported by the Natural Science and Engineering Re- 
search Council of Canada, and the research of CSL is also supported by the 
Quebec Department of Education. This investigation was initiated during a 
visit of CSL to TRIUMF. He wishes to thank the Theory Group there for its 
hospitality. He also wishes to thank Keith Dienes for informative discussions. 
JNN would also like to thank Professor T.K.Lee of the National Center for 
Theoretical Science, Taiwan for his hospitality where part of this research is 
done. 



References 

[1] T. Kaluza, Math. Phys. Kl. (1921) 966; O. Klein, Phys. Z. 37 (1926) 895. 
[2] J. Polchinski, Phys. Rev. Lett. 75 (1995) 4724. 

[3] N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B429 (1998) 
263; Phys. Rev. D59 (1999) 086004; I. Antoniadis, N. Arkani-Hamed, S. 
Dimopoulos, and G. Dvah, Phys. Lett. B436 (1998) 257; L. Randall and 
R. Sundrum, Phys. Rev. Lett. 83(1999) 3370, 4690. 

[4] J.C. Long, A.B. Churnside, and J.C. Price; CD. Hoyle et. al., 
Phys. Rev. Lett. 86 (2001) 1418. 

[5] See, e.g., J.M. Cline, C. Grojean and G. Servant, Phys. Rev. Lett. 
83(1999) 4245; P. Binetruy, C. Deffayet, U. EUwanger and D. Langlois, 
Phys. Lett. B477 (2000) 285; C. Csaki, M. Graesser, L. Randall and J. 
Terning, Phys. Rev. D62 (2000) 045015, and references therein. 

[6] K.R. Dienes, E. Dudas, and T. Gherghetta, Nucl. Phys. B557 (1999) 
25; K.R. Dienes and I. Sarcevic, Phys. Lett. B500 (2001) 133-141; D.O. 
Caldwell, R.N. Mohapatra, and S.J. Yellin, |hep-ph/001Q353| ; 0102279. 



15 



[7] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, and J. March-Russell, |hep- 



ph/9811448t G. Dvali and A.Yu. Smirnov, Nucl. Phys. B563 (1999) 63; 
R.N. Mohapatra, S. Nandi, and A. Perez-Lorenzana, Phys. Lett. B466 
(1999) 115; Nucl. Phy. B593 (2001); R. Barbieri, P. Creminelli, and A. 
Strumia, Nucl. Phys. B585 (2000) 28; R.N. Mohapatra and A. Perez- 
Lorenzana, Nucl. Phys.B 576 (2000) 466; Nucl. Phys. B593 (2001) 451; 
N. Cosme, J.-M. Frere, Y. Gouverneur, F.-S.Ling, D. Monderen, and V. 
Van Elewyck, Phys. Rev. D63 (2001) 113018. 

[8] A.E. Faraggi and M. Pospelov, Phys. Lett. B458 (1999) 237; A. Das 
and O.C.W. Kong, Phys. Lett. B470 (1999) 149; G.C. McLaughlin and 
J.N. Ng, Phys. Lett. B470 (1999) 157; B493 (2000) 88; A. Lukas, P. 
Ramond, A. Romanino, and G.G. Ross, Phys. Lett. B495 (2000) 136; 
hep-ph/0011295 . 



[9] J.N. Bahcall, P.I. Krastev, and A.Yu. Smirnov, |hep-ph/0103T79 



[10] J. Bonn et al, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 273. 
[11] H.V. Klapdor-Kleingrothaus, |hep-ph/01023T9| . 

[12] The Super-Kamiokande Collaboration, Phys. Rev. Lett. 85 (2000) 3999- 
4003 (|hep-ex/0009001| ); H. Sobel, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 
127; Y. Suzuki, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 29. 

[13] G.I. Fogli, E. Lisi, and A. Marrone, Phys. Rev. D63 (2001) 053008. 

[14] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theo. Phys. 28 (1962) 870. 



16 




Figure 1: Graphical solutions of the characteristic equation (H). The left hand 
side of (1^) is represented by the solid curve; the right hand side is represented 
by the dashed curves, with (7711,7712,^3) = (—1.8,-0.2,0.3). The one with 
the short dash has the larger coupling constant d than the one with the long 
dash. Note that the solution A moves up and right along the solid curve for 
increasing d when the value on both sides of (|^) is positive, and moves down 
and left along the solid curve when the value is negative. 
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